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Abstract —A recent claim by Lieu et al that beam splitter 
intensity subtraction (or homodyne with one vacuum port) 
followed by high resolution sampling can lead to detection of 
brightness of thermal light at the shot noise limit is reexamined 
here. We confirm the calculation of Zmuidzinas that the claim 
of Lieu et al was falsified by an incorrect assumption about the 
correlations in thermal noise. 

I. Introduction 

Recently, (I] proposed a method of improving the sen¬ 
sitivity of radio telescopes, based on using a 50:50 beam 
splitter and measuring the difference signal between the two 
output beam intensities. The original motivation of III was 
based upon the supposition that the intensity difference has 
shot noise fluctuations given by a simple Poisson process 
with its mean subtracted away, and a variance equal to the 
mean photon rate of the(presumed stationary) incident beam, 
although the statistics of the noise distribution presented there 
were derived from the quantum theory of chaotic light and 
the higher moments were found to be slightly different from 
Poisson. Subsequently, it was pointed out by |2l that the 
higher moments of were erroneous, due to an invalid 
assumption about the absence of correlations between non¬ 
overlapping time intervals. Below, we present an improved 
and simplified version of the calculation of H] that reaches the 
same conclusion as 11, viz. although it is possible to use the 
split-beam technique to achieve essentially the same accuracy 
as that of a direct measurement of the incoming signal, it is 
much harder (if at all possible) to do significantly better. 

II. An Improved Calculation 

This is firstly an alternative calculation to the ones in IT| 
and 0) of the degree of accuracy attainable by a difference 
measurement on a split beam, and secondly a discussion of 
the effect in frequency space via a discrete Fourier transform. 
The results are essentially the same as those of Zmuidzinas, 
i.e. the split-beam technique achieves more or less the same 
accuracy that one can get from a direct measurement of the 
brightness of the original beam. 

Ahead of the formal treatment, it may be useful to seek a 
heuristic understanding of the difference between HI and 0. 
If light comprises only shot noise, the fluctuations in direct 
and homodyne measurements will in principle both look the 
same. But, as is usually the case for chaotic light, there are 
classical phase noise fluctuations as well, and the shot noise 
will then exhibit a time dependent mean and variance as its 


amplitude varies together with the classical intensity noise 
in tandem. According to 0 and the calculation here, this 
correlation between the shot noise variance and the classical 
noise intensity is not expected to be removed by the beam 
splitter, as illustrated in Figure 1. 0, on the other hand, 
asserted that the appearance of the homodyne difference signal 
remains like simple shot noise, i.e. either there was no such 
con'elation in the incident beam to begin with, or the beam¬ 
splitter removed the effect. 


A. Direct measurements with incident beam 

So long as we are dealing with a narrow bandwidth, it is 
more convenient to work with the Fourier transforms of the 
annihilation and creation operators, 

d{t) = —^= f duj d\t) = -^= f dw 

V 27r j s/2'k J 

( 1 ) 

They satisfy the commutation relations 

[d{t),d{t')]= 5{t-t'). (2) 


For a chaotic beam with Gaussian frequency profile, centted 
on Wo and with bandwidth 1/t, we have 

{a'{t)d{t')) = nof{t - t'), (3) 

where f(t) is given by 

f{t) = (4) 

T 

Moreover, the intensity is simply woa^(f)a((:), but in the 
narrow-band case, it is simpler to remove the factor of wq, 
and talk instead about 


J{t) = d\t)d{t), (5) 


which represents the number of photons arriving per unit time. 
It follows immediately that 

{J{t)) = -• (6) 

T 

Next, we examine the covariance function 

cov( J(f), J(f')) = (Jit) Jit')) - (J)". (7) 


Now 

{Jit)Jit')) = {d\t)d{t)a'it')dit')) 

= {a\t)dit)){d\t')ait')) + 
{d\t)ait')){a{t)d\t')). 


( 8 ) 





j 



Fig. 1. Theoretically expected noise characteristics of stationary thermal 
radiation as viewed directly and via a beam-splitter. Note the squared intensity 
difference still exhibits classical bunching noise that’s con'elated with the 
direct intensity time series Sj. Additionally Dj is also expected to have more 
shot noise. 


The first term here clearly cancels the last term in Q, so we 
find 


cov{J{t),J{t')) = nl\f{t - + —S{t - t'), (9) 

T 

where the final term comes from the commutator [a(f), a^{t')]. 
If we define the average flux over a short time interval as 


Jrit) 




dt' 


t-T 


( 10 ) 


then we find 


where 


va,T{JT{t)) 



F 





( 11 ) 


( 12 ) 


Note that for T t, we may replace / in the integrand 
by 1/r, so F{T/t) « T/t. The relative uncertainty in the 


measurement of Jj- is given by 

var(Jr(f)) ^ 1 L V J_1 

(JT(f))2 t[ \t) noj ’ 


(13) 


or 


var(Jx(^)) 


IH-for T<t. 

noT 


(14) 


On the other hand, if we measure for a much longer time 
T = NT, we must use the limiting value of f{x) for a; 3> 1, 
namely yTr. So we have 


var( Jr(f)) 


T:>t. (15) 


B. Dijference signal for split beam 

In a 50:50 beam splitter, it is useful to consider a second 
input beam, which is in fact in its vacuum state. Let us 
represent the annihilation and creation operators of that second 
input by b{f),F{t). Then for the two output beams we have 
annihilation operators 


c = 


1 


(a + ib), 


d 




(16) 


Note that c and d each satisfy the commutation relations (|^, 
together with ([^ but with no replaced by no/2. Moreover, 
[c, Jl] = 0. 

One might perhaps worry that using h{f) rather than b{u:), 
with the replacement of factors of oj by ujq, which is justified 
for the narrow-bandwidth case, might be inadmissible for the 
vacuum contribution. However, if one retains the factors of w, 
they will be converted to time derivatives that will ultimately 
act on other factors that are limited by bandwidth, and the 
leading contributions will be given quite accurately by the 
replacement of w by ojq, so this is probably not a serious 
problem. 

The quantity we are particularly interested in is the differ¬ 
ence signal, the difference between the numbers of photons 
arriving in the two output channels. This is given by 


D{t) = c\t)c{t) — S{t)d{t). (17) 


Substituting from we see that this quantity may be written 
D{t) = ia^{t)b{t) — iP (t)d{f). (18) 


Obviously, its expectation value is zero: 

{D{t)) = 0. (19) 

The factorization between d and b operators makes this a very 
convenient form to use. For example, in computing the two- 
time function, we see that 

( 20 ) 

and because the b input is in its vacuum state, the second term 
vanishes, while in the first, (S(f)6^(f')) = S{t — t'). Thus we 
find 

cov[D{t),D{t')) = {D{t)D{t')) = —5{t - t'). (21) 

r 

























So the measurement of the variance of D provides a way of 
measuring np. 

Of course, any measurement will take up a finite time 
interval. We suppose that the total available time T is divided 
up into N small segments of duration T, and define the average 
flux in the jth interval as 


Dj = - 

j j. 


rjT 


dtD{t), 


Putting these expressions together and substituting into (25 1 , 
we find 


(5(fi2)^(f34)n-o ( ^ + l/(^13)P 


( 22 ) 


'0-l)T 

where we assume T ^ r, so that 

variDj) = cov{Dj,Dk) = 0, {j ^ k). (23) 

Now to estimate the accuracy of the measurement we can 
make, we need to compute the expectation value 
However, for use later we consider the more general case 

1 

(DjDkDiDm) = dti dt2 

J Ju-1)T J(k-1)T 

plT pmT 

/ dt^ / dt4 

{b{h)b{t2)D{h)b{u)). (24) 

When we substitute from (0, each term in the expectation 
value can be written as a product of an expectation value of d 
and operators, and one of b and operators. Moreover, the 
latter vanish if they have a 6 on the right or a b on the left, 
and there must be equal numbers of each of the two terms in 
(18 1 containing b and b operators. So there are just two terms 
remaining: 

{b{ti)b{t2)b{tb)b{t4)) = {b {tx)b {t2)d{t'i)d{t4)) 

{h{tx)b{t2)b\tb)V' (tb) + 

(a’^' {ti)d(t2)b bbaiU)) 
{b(ti)b {t2)b{t3)b (tb) ■ (25) 

Now, with the abbreviation tjk = tj — tk, 

(b{ti)b(t2)b{t3)bbb) = S{ti3)S{t2b + S{tibS{t23), (26) 

while 

{b{ti)b {t2)b{t3)b bb) = ^(fl2)<J(f34), (27) 

so clearly the result will only be nonzero when the indices 
(j, fc, I, m) are equal in pairs. 

We also note that 


{Dbi)Db2)DbbDbb) = 
1 
T' 

1 
T- 
1 


<5(^r3)^(^24)?^o I 2 


(5(^14)^(^23)?^i 


l/(tl2)P 


Sbl2)^b23)Sb3b 


no 


(30) 


Note the symmetry of this expression under permutations of 
{1,2, 3,4}, which results from the fact that the different bj 
operators commute with each other. 

Then, integrating over short time intervals, and assuming 
that r <C T, we find 


^b jbkbib'b) — djkSir, 


Th 


T2t2 


2(1 + 6 "*?'/^") + 


^jmSklrT^b-2 (1 + e ) + 




no 


djkdkblm rr^n ; 

1 


(31) 


where tjk = (j — k)T. 

Now, to find the covariance of Dj and Df, we set k = j 
and m = I, and remove the first of the seven terms in ( (3T] l, 
which is cancelled by the product of expectation values. This 
yields 




'ji 


T2t2 


Ojl 


no 

T^t 


(32) 


The first term alone gives the covariance when j b For 
j = I we find 




(33) 


(34) 


Thus the fractional etTor is given by 

var(Zl2) ^ .r 

{D]b noT 

This is comparable with ( [l4| ) but larger (when uqT ^ r) by 
a factor of 5. 

Of course, as before we can do better by observing for a 
longer time. In particular, we can form the sample mean 


{b bi)b b2)db3)bb)) = bo[fbi3)fb2b + fbibf(t23)]- 

(28) 

while 

{b bi)db2)b b3)db4:)) = nl[fbi2)fb3b + 

fbi‘i)fb32)] + 
?^o/(^14)(5(^23)■ (29) 


1 


N 


D^ = bynl 

N ^ ^ 


(35) 


Clearly, 


i=i 


N 


var(D2) = ^ ^ cov{D],Df). 


3,1 = -^ 


(36) 











When we substitute from ( |3^ , in the first term, we can convert 
the sum over j — I to a Gaussian integral: 


j 

Thus we obtain 

var(i:>^) ^ T \ 

{D^)2 N \ T ^ ^ noTj ' 


Tj T 


(37) 


(38) 


The dominant term here, when T t uqT is the first. 
This reproduces precisely the dominant term in the direct 
measurement error, ( [T5] l. So in this case we can do as well 
as the direct measurement, but unfortunately no better. 


C. Frequency domain measurements 

There is, however, another way of dealing with the infor¬ 
mation, in terms of the discrete finite Fourier transform of the 
signal. Let us define 


1 


N 




^ — 27 ripj /N 


(p = 0,...,7V-l). (39) 


j=t 


The frequency corresponding to Kp is ujp = 2'Kp/NT. Note 
that Kp is not real; in fact K* — K-p = K^-p- Obviously, 
{Kp) = 0, and, from (231, 


{KpkD = 


riQ 

Tr’ 


{kpKl) =0, (p ^ q), (40) 


very similar to the expressions for the variance and covariance 
of the Dj. Thus for each p the value of provides an 

estimate of uq. 

Next, we look at the uncertainty of these estimates. To do 
so, we need the covariance of and |iTqp. We start from 
and apply a discrete Fourier transform to each of the 
variables. This yields 

9 

-N ^ /s ^ 

{KpKgKrKs) — ^p-\-q+r-\-s[/ji2^2 Gp-\-q Sp-\-r 

Gp+r + 6p+s + Gp+s) + (41) 

where (5„ = 1 if p = 0 mod N, otherwise 0, and 


Gp = — ^ /r^ 


s/^T _ 

~NT^ 


{pttt/NT)^ 


(42) 

To find the covariance of |iTpP and we have to set 

q = —p, s = —r, and again subtract the first term, to obtain 


cov(|iTpp, \Kr 


''0 


+ Gp-r) + 


{Sp+r + 5p-r + Go + Gp+r 


no 


p-^' ' NTh 

In particular, setting p = r we have 


var(|iTpp) = 


'•0 


T2t2 


NT 


1 + ^p + 2 


„-(2p7rr/Afr)' 


NT 


(43) 


no 

NT^t 


(44) 


The fractional uncertainty measure is now 

var(|iTp|2) ^ ^ ^ ^ , V^r 

{kpkl)^ + + + ^noNT- 

(45) 

This already shows that when T t NT, the variance of 
each |iTpp is less than of each Dj. 

As with the Dj, we can form the sample mean of these 
measurements, 

_ 1 ^-1 

=N^ 

p—0 

so that 

_ 1 N-i 

var(|A:|2) = — covdATplMATpn. (47) 

p^r—0 


For the terms involving Gp we can again convert the sum to 
an integral to give 




Thus 


var(|Ar|2) 


^0 ( 4 +^]+^}^ 

NT^t^ \tJ NT^K 


whence 


var(|iT|2) 

{kpkl)^ 


1 

N 




(48) 


(49) 


(50) 


Note that this is precisely the same as the result for D"^. 
This might seem surprising, given that the variance of each 
individual \K‘^\ is less than that of each i9|. However, it arises 
from the fact that there is some correlation between all the 
|iTp|, whereas for the it is limited to time differences less 
than about r. The situation is depicted in figure (|II-C 1 . 



Fig. 2. The Fourier power spectrum of the time series Dj has both random 
and systematic uncertainties, with the latter in the form of a normalization 
error that shifts the entire (flat) spectrum up or down as measurements of the 
same thermal source are repeated. The origin of this systematic effect, which 
for chaotic thermal light dominates the random noise when the spectral data 
are averaged over the full range of available frequencies, is photon bunching 
noise (also known as classical noise. 








































D. Fourier transforming 

We could also envisage a different use of the finite discrete 
Fourier transform. Let us define the quantity Ej = D^, which 
of course satisfies 

( 4 ) = ( 51 ) 

Then we may define its discrete Fourier transform, say 
1 ^ 

^ E {p = 0,...,N-l). (52) 

It follows at once that 

(4) = VN^Sp. (53) 

So in this case none of the Fourier components, except the DC 
(which apart from a normalization factor is the same thing as 
D^), can provide a measure of uq. So this does not seem a 
very profitable avenue of enquiry. 

III. Conclusion 

We revisited a recent claim IT] that beam splitter intensity 
subtraction (or homodyne with one vacuum port) followed by 
high resolution sampling can lead to detection of brightness 
of thermal light at the shot noise limit. We are able to confirm 
the result of la that the claim of |l] was falsified by their 
incorrect assumption about thermal noise. 

We should emphasize, however, that there has so far not 
been any experimental comparison of direct against homodyne 
brightness measurement of the same thermal state. Thus the 
quantum field theoretic predictions remain to be verified in the 
laboratory. 
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